Knudsen-to-Hydrodynamic Crossover in Liquid 3 He in a High-Porosity Aerogel 
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We present a combined experimental and theoretical study of the drag force acting on a high 
porosity aerogel immersed in liquid 3 He and its effect on sound propagation. The drag force is 
characterized by the Knudsen number, which is defined as the ratio of the quasiparticle mean free 
path to the radius of an aerogel strand. Evidence of the Knudsen-hydrodynamic crossover is clearly 
demonstrated by a drastic change in the temperature dependence of ultrasound attenuation in 98 % 
porosity aerogel. Our theoretical analysis shows that the frictional sound damping caused by the 
drag force is governed by distinct laws in the two regimes, providing excellent agreement with the 
experimental observation. 

PACS numbers: 67.30.eh, 47.56.+r, 47.61.-k 



The understanding of complex flow phenomena, such 
as in porous media, is a challenging goal of fluid physics. 
The pioneering work by Darcy [l( established empirically 
a linear relation between the fluid flux of laminar viscous 
flow in porous media and the total pressure drop in bal- 
ance with viscous drag force. The Darcy law provides 
a constitutive relation in hydrodynamics and has been 
applied to various kinds of flow phenomena such as wa- 
ter flow in aquifers, blood flow in vessels, and petroleum 
flow through sandstone and gravel ■ In general, how- 
ever, the hydrodynamic description breaks down when 
the mean free path (mfp) of fluid particles exceeds sub- 
stantially the characteristic dimension of the porous ma- 
terial. In this case, the fluid mass is carried by the 
Knudsen diffusion [5], i.e. , successive ballistic flights of 
individual particles followed by scattering at walls. The 
Knudsen diffusion, as well as the viscous flow, is of partic- 
ular interest for gas transport in nanofabricated systems 
0] > which have recently attracted significant attention 
in the context of nanofluidics Q. The investigation of 
the mfp dependent flow regimes is also of fundamental 
importance towards a unified understanding of the trans- 
port properties in complex porous materials existing in 
nature. 

Due to the unique porous structure of aerogel and the 
fascinating low temperature properties of liquid 3 He, the 
3 He-aerogel system presents an opportunity that allows 
investigations of a wide range of physical phenomena. 
The superfluid phase of 3 He in highly porous aerogel is 
an analog of unconventional superconductors with pair- 
breaking impurities and has been the subject of 
intense study since the first observations 0, [l(| of su- 
perfluid transitions in this system. In the normal liquid, 



one can study the Fermi liquid behavior ranging from 
transport properties ll|, Il5l4l9| to diverse flow phenom- 
ena [201 ] in the presence of impurity-like disorder with 
a topologically antithetical geometry to the other con- 
ventional porous media. In this Letter, we present the 
first experimental and theoretical study on the Knudsen- 
hydrodynamic crossover in the normal liquid 3 He con- 
fined in 98% porosity silica aerogel. 

Liquid 3 He is a paradigm of Fermi liquids in which 
the transport property is governed by low energy ther- 
mal excitations called quasiparticles. Since the collision 
processes between the quasiparticles are restricted by the 
Pauli exclusion principle, the quasiparticle mfp, I, can be 
much longer than the interparticle distance. The mfp in- 
creases with decreasing temperature T, following a 1 /In- 
dependence [21| . and approaches ~ 10 /jm near the su- 
perfluid transition temperature, T c ~ 2 mK (see Fig.QJ. 
On the other hand, aerogel, an interconnected network of 
thin strands of aggregated silica particles, introduces ad- 
ditional length scales: the strand diameter 2a, the inter- 
strand distance and the geometrical mean free path 
l g . Specifically, in 98% porosity aerogel, 2a ~ 3 nm, 
C « 25 - 40 nm, and finally l g « 100 - 150 nm. There 
is no well-defined pore structure in aerogel in contrast to 
the nanochannels used in the previous works 0, 0] on the 
Knudsen-hydrodynamic crossover. 

To understand the effect of the high porosity aerogel 
on the mass transport in liquid 3 He, we performed exten- 
sive ultrasound transmission measurements over a wide 
temperature range up to 200 mK. Based on the exper- 
imental results, we shall show that the transport prop- 
erties in aerogel are characterized by two length scales 
of aerogel, the strand radius a and the geometrical mfp 
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FIG. 1: Characteristic length scales of 98% porosity aerogel 
and the temperature dependence of the quasiparticle mean 
free path for three pressures, I oc 1/T 2 in pure liquid 3 He. The 
mfp was extracted from the experimental data of viscosity 
[221 ] . The whole temperature range is divided into three trans- 
port regimes classified by the hydrodynamic, quasi-Knudsen, 
and full Knudsen regime. For T < 15 mK (full Knudsen 
regime) the quasiparticle mfp is limited by the geometrical 
mfp l 3 . Above ~ 15 mK, the mutual quasiparticle collision 
becomes dominant, but the drag force is still in the Knud- 
sen regime as far as I > a (quasi-Kundusen regime). Around 
150 mK (hydrodynamic regime), where I ~ a, the Knudsen- 
hydrodynamic crossover of the drag force occurs. The figures 
on the right show a schematic illustration of two extreme flow 
regimes. 



l g , separating three transport regimes as summarized in 
Fig.ffl 

The acoustic cavity used in this experiment is formed 
by two longitudinal LiNbC>3 transducers used as a trans- 
mitter and a receiver. The transducers are separated 
by a Macor spacer maintaining a 3.05 ± 0.02 mm par- 
allel gap. The pair of transducers were carefully se- 
lected to achieve the best match in the fundamental fre- 
quency (9.5 MHz) and the frequency response. The 98% 
porosity aerogel sample was grown between the trans- 
ducers in situ to prevent any irregular reflections at the 
transducer- aerogel boundaries. Temperature was deter- 
mined by a melting pressure thermometer for T > 1 mK 
and a Pt NMR thermometer for T < 1 mK. The pres- 
sure of the cell was monitored throughout the measure- 
ments by a capacitive pressure sensor incorporated as a 
part of the experimental cell. A commercial spectrom- 
eter, LIBRA/NMRKIT II (Tecmag Inc., Houston, TX), 
was used to excite the transmitter and observe the prop- 
agated sound signal from the receiver. Alfis pulse was 
sent to the transmitter and the signal from the receiver 
was amplified with a preamplifier (Miteq AU1114) be- 
fore being transferred to the spectrometer. A detailed 
description of the experimental setup can be found else- 
where [l5j |. 

In Fig. [51 the receiver signal outputs acquired at 29 bar 
are plotted for selected temperatures. The signals for all 




FIG. 2: Receiver signal outputs at 29 bar for selected tem- 
peratures. The superfiuid transition in aerogel is demarcated 
by an arrow. 



temperatures rise around 8 fxs. Below 0.5 mK, in the low 
attenuation regime, one can see the three echoes following 
the initial peak. From the time-of-flight measurements, 
the sound velocity c was obtained to be c = 330 m/s 
at 29 bar. The temperature dependence of the sound 
velocity is less than our measurement error for the whole 
temperature range presented in this work. 

The relative sound attenuation, a, was determined us- 
ing two different methods. One was to use the peak value 
of the initial signal and the other was to use the inte- 
grated signal of the initial peak from the rising edge to 
23 /its. Both methods produced consistent results within 
1% in the low attenuation regime and within 5 % in 
the high attenuation regime. The relative attenuation 
was then converted into the absolute attenuation using 
a reference value at 0.4 mK where 3 echoes were clearly 
observed owing to the low attenuation in the superuid 
phase. The peaks of the echoes exhibit an exponential 
decay, enabling us to determine accurately the absolute 
attenuation. In Fig. [3j the attenuation data obtained 
by the second method are plotted as a function of tem- 
perature T. The data at three different pressures are 
shown above the superfiuid transition in aerogel. The 
attenuation a is almost independent of T near T c but 
decreases gradually with increasing T. This behavior is 
in good agreement with those observed in previous ul- 
trasound attenuation measurements However, our 
higher temperature measurements revealed a new behav- 
ior which exhibits a minimum around 50 mK followed by 
a rather rapid increase in attenuation with temperature. 

As the first step of analysis of the experimental re- 
sults, let us estimate l g in 98 % aerogel. The simplest 
estimate can be made by modeling the aerogel as a di- 
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FIG. 3: Ultrasound attenuation (uj/2n = 9.5 MHz) in the 
normal fluid above the aerogel superfluid transition. The solid 
symbols represent the experimental data at three different 
pressures. The solid lines are the results of fit to the data 
using Eqs. (Q-@ with the Cunningham constant A = 6/5. 
The fitting parameters are (lti,la) = (115.4, 157.2) nm for 10 
bar, (108.7, 147.1) nm for 21 bar, and (100.0, 135.4) nm for 29 
bar. The dashed lines are the sound attenuation calculated in 
the Knudsen (Qoo) and hydrodynamic (qq) limit (see text). 



the geometry-limited value of (l/5)npF l a - Similar satu- 
ration behavior is observed for the spin-diffusion coeffi- 
cient of 3 He in 98 % aerogel [Uj]. 

Now we discuss the sound propagation in aerogel. The 
sound frequency uj = 2n x 9.5 MHz used in our mea- 
surements should be compared to the effective quasi- 
particle scattering rate, l/r c fj = VF(1/Z a + 1/0- The 
Fermi velocity Vf in liquid 3 He is ~ 50 m/s [25j. Since 
1/T c ff > vf/Io, > 10 s s -1 , our ultrasound measurements 
are in the collision-dominated regime, uj < l/r fj. There- 
fore, the sound propagation can then be described by the 
Navier-Stokes equation 

pd t v = -VP + r lcS {Av + ^\7divv)-F, (2) 

where p = ran is the 3 He mass density, v the velocity 
field of liquid 3 He, P pressure, and F the drag force (per 
unit volume) exerted on aerogel. This drag force gives 
rise to the coupled oscillations of liquid 3 He and aero- 
gel 0, [It}- Since the sound wavelength of interest is 
much longer than the distance between aerogel strands, 
the aerogel oscillations can be described within a contin- 
uum approximation, 



PadtU a = p a C^Au a + F, 



(3) 



lute random distribution of spheres with radius a (sphere 
model). Noting that the sphere density n a is related to 
the aerogel porosity </> by (4/3)7ra 3 n a = 1 — 0, we ob- 



tain L 



1/Tra 2 n a = (4/3)o/(l 



Substitution of 



a = 1.5 nm and <f> = 0.98 yields l g = 100 nm. 
realistic model of randomly oriented cylinders 
leads to a similar value l g = 2a/(l — <p 



A more 
231 ] also 
150 nm. These 
estimations are consistent with those from an empirical 
formula extracted from computer simulations based on a 
diffusion limited cluster aggregation model 24[. In Fig. 
[U l g is compared with I oc 1/T 2 in pure liquid 3 He. We 
observe l g ~ I at around 15 mK. 

At low temperatures where I > l g , the transport prop- 
erties are governed by the geometry of aerogel struc- 
ture. This point is clearly demonstrated by regarding 
the quasiparticle-quasiparticle and quasiparticle- aerogel 
collisions as two independent scattering sources. Then, 
for example, the effective viscosity coefficient in the 3 He- 
aerogel system has the form 



VeS = \npF{l/la + 1/0 1 = 7 



np F la 

5 1 + tvt; 2 



(1) 



with n being the 3 He number density, pp the Fermi mo- 
mentum, l a a temperature independent mfp of the or- 
der of l g , and T a the temperature at which I = l a . At 
high temperatures, T 3> T a , Eq. (QJ is reduced to the 
conventional gas-kinetic formula for a pure Fermi liq- 
uid. As temperature decreases, the quasiparticle scat- 
tering processes become gradually dominated by the 
aerogel-quasiparticle collisions, so that 7y e ff saturates to 



where u a , p a — 0.044 g/cm 3 , and c a ~ 50 m/s are the 
displacement vector, the mass density, and the speed of 
longitudinal sound in aerogel itself, respectively. 

To proceed further, we need the constitutive relation 
for F. Let Kn be the Knudsen number defined by Kn = 
I /a oc 1/T 2 . Since the smallest length scale of aerogel is 
the strand radius a, the hydrodynamic description of the 
drag force is justified for Kn < 1. In the hydrodynamic 
limit, Kn — > 0, the drag force density obeys the Darcy 
law 



F = -v rel . 

K 



(4) 



Here v rc \ = v — dtu a is the relative velocity field of liquid 
3 He and aerogel, r\ = (l/5)nppl the viscosity of pure 
liquid 3 He, and k the permeability. In the sphere model 
for aerogel, the Darcy law is equivalent to the Stokes law 
Fo/n a = 6Ttar]v lc i, and then the permeability takes the 
form, k — al g /6. 

In the opposite limit, Kn — > oo, where the Darcy law 
breaks down, the drag force is determined by the colli- 
sions between aerogel and individual quasiparticles (see 
Fig- ID- Ichikawa et al. [l6[ discussed the drag force on 
aerogel by modeling it as randomly distributed impuri- 
ties. By calculating the momentum transfer from liquid 
3 He to aerogel using the impurity-quasiparticle collision 
integral, they obtained the following expression for the 
drag force density: 



npF 



(5) 
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where Z tr is the transport mfp. Equation ([5]) has an an- 
ticipated form of the drag force density in the Knudsen 
limit, Kn — > oo. Specifically, consider a situation in which 
an object with a cross section a is moving with a constant 
velocity — v Te \ in "gas" of quasiparticles with momentum 
pp. The body would experience collisions with the quasi- 
particles of number ncr|i> r ei| in unit time and experiences 
a force ~ pFnav T0 \. Dividing it by the total volume V 
and replacing a/V by 1/Z tr , we reach Eq. ((5]). It should 
be noted, however, that since /t r depends on the details 
of the object, e.g. , its shape and the irregularity of the 
surface, it is difficult to estimate Z t r accurately for com- 
plex porous materials like aerogel. We should therefore 
regard Ztr as a phenomenological parameter comparable 
to the geometrical mfp l g . 

In order to obtain the expression for the drag force 
for arbitrary values of Kn, we put F = Fq/(1 + AKn) 
and determine the dimensionless constant A so that the 
Knudsen-limit formula Eq. ([5]) is recovered. A similar 



procedure was used by Cunningham |26| in his theory 
for the drag force on a sphere immersed in classical flu- 
ids, and the resulting Stokes-Cunningham formula was 
applied successfully to Millikan's oil-drop experiment in 
air [27( 1 . In the present problem, the Cunningham con- 
stant A takes the form of A — al tr /hn. For the sphere 
model, we have A = 6Zt r /5Z g , which illustrates that A is 
a constant of order unity. For the analysis of the experi- 
mental results, it is convenient to write F in the following 
form: 



1 



1 



np F 



1 + I /AKn 



l + T 2 /AT 2 l t 



-v Te \. 



(6) 



Here, T\ is the temperature at which Kn = 1. The 
crossover temperature T\ is well above T a , as can be seen 
in Fig. [Q 

The Cunningham-type formula implicitly assumes a 
smooth transition of the drag force between the hydro- 
dynamic and Knudsen regimes. Virtanen and Thuneberg 



29j recently reported a detailed numerical study of 
the drag force on a vibrating wire immersed in 3 He- 4 He 
mixtures. Although the system is different from that of 
present interest, a common feature can be found in their 
results, showing a smooth transition in the dissipativc 
part of the mechanical impedance Z at low frequencies. 

Having obtained the system of equations, we can now 
analyze the experiment. Considering the longitudinal 
sound of the plane wave form oc e l<yq ' r ~ ult ^ and taking ac- 
count of an experimental fact that \c a q/oj\ 2 ~ 0.02 <^ 1, 
one obtains 



where c = cj/Req = ci/(l + p a / p) 1 ^ 2 and c\ are the 
sound velocities in the presence and absence of aero- 
gel, respectively. The sound velocity c is reduced from 
ci by a factor (1 + pf p a )~ x l 2 due to the extra inertia 
from the dragged aerogel. Using ci = 398 m/s and 
p = 0.114 g/cm 3 at 29 bar [Hj], we obtain c = 338 m/s, 
in good agreement with the experiment. c\ and there- 
fore c are almost independent of temperature, which also 
agrees with our observation. In Eq. J?} for sound atten- 
uation, a v and a / are associated with the viscous damp- 
ing and the frictional damping caused by the drag force, 
respectively. The theoretical calculation represented by 
the solid lines in Fig. [3] shows excellent agreement with 
the experiment using the fitting parameters, In and l a as 
listed in the figure caption and A = 6/5. The role of the 
Cunningham constant A is to shift the crossover temper- 
ature in Eq. ([9]). The best fit to the experimental data 
is achieved by taking A to be around 6/5, corresponding 
to the sphere model with Z tr = l g . 

To gain insights into the observed behavior, we plot in 
Fig. [3] the sound attenuation in the two limiting cases for 
29 bar: the Knudsen (aioo) and the hydrodynamic (ao) 
regimes by taking the limits of T\ — > 00 and T\ — > 0, 
respectively in Eq. (9). The difference between the two 
limits comes mainly from the different T-dependence of 
the frictional damping aj. In the full Knudsen regime 
where T T a and effectively T\ — > 00, both a v and 
a/ are temperature independent. This simply reflects 
the fact that the effective quasiparticle mfp is limited by 
the geometrical mfp l g . The decrease of the Knudsen at- 
tenuation above T a is due to the decrease of a v . The 
low temperature experimental data are well described 
by the Knudsen behavior. However, as T increases, a 
significant deviation occurs especially around T = T\. 
With increasing T further above Ti, the experimental 
data approach asymptotically the hydrodynamic behav- 
ior, a w ao a T 2 . The crossover temperature T\ in- 
creases with decreasing pressure (see Fig. [T]) due to the 
pressure dependent quasiparticle mean free path. Con- 
sequently, this effect shifts the attenuation minimum to 
higher temperature for lower pressure as observed in the 
experiment. 
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